Quantum Stress Focusing in Descriptive Chemistry 
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We show that several important concepts of descriptive chemistry, such as atomic shells, bonding 
electron pairs and lone electron pairs, may be described in terms of quantum stress focusing, i.e. 
the spontaneous formation of high-pressure regions in an electron gas. This description subsumes 
previous mathematical constructions, such as the Laplacian of the density and the electron local- 
ization function, and provides a new tool for visualizing chemical structure. We also show that the 
full stress tensor, defined as the derivative of the energy with respect to a local deformation, can be 
easily calculated from density functional theory. 
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Atomic shell structure, electron pair domains, tt- 
electron subsystems, etc., are common concepts in de- 
scriptive chemistry and play a significant role in modern 
electronic structure theory. These concepts help us to 
visualize the bonding between atoms in terms of small 
groups of localized electrons (e.g. two electrons of oppo- 
site spin in a simple covalent bond) and therefore play 
an important role in predicting new molecular structures 
and in describing structural changes due to chemical re- 
actions. 

A precise quantitative description of small groups of 
localized electrons has been sought for a long time [H, [2| , 
but none of the solutions proposed so far is completely 
satisfactory. The natural candidate - the electronic den- 
sity n(r) - has a clear physical meaning, but fails to re- 
veal quantum mechanical features such as atomic shell 
structure ^ or the localization of electron pairs of oppo- 
site spin in a covalent bond. The Laplacian V^n of the 
density [H, Q provides a better way to visualize molecu- 
lar geometry, but lacks a clear physical significance and 
fails to reveal the atomic shell structure of heavy 
atoms. Recently, a very useful indicator of electron lo- 
calization has been constructed Q from the curvature 
of the (spherically averaged) conditional pair probability 
function Paaiic, r') (the probability of finding an electron 
of spin cr at r given that there is another electron of the 
same spin orientation at r') evaluated at r = r'. In the 
simplest approximation this curvature is proportional to 



(1) 



where T„{r) = i|V'!/'/(T(r)p is the non-interacting ki- 
netic energy density of cr-spin electrons (atomic units 
m = h — — I are used throughout), V'io-(r) are the 
occupied Kohn-Sham orbitals of spin a, and rio- is the 
(T-spin electron density (n = + ni). From Da one 



constructs the electron localization function (ELF)[^ 
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?7a(r) = 



l + (i?,/TjF)2' 



(2) 



where rj^ — {3/10){6T:^)^^^na^^ is the Thomas-Fermi 
kinetic energy density of tr-spin electrons. The ELF pro- 
vides excellent visualization of atomic shells as well as 
a quantitative description of valence shell electron pair 
repulsion theory of bonding and it has been recently 
extended to excited and time-dependent states 0] • Nev- 
ertheless, it has two defects: (1) it remains a mathemat- 
ical construction of dubious physical significance Q and 
(2) it is difficult to calculate ELF beyond Q the lowest- 
order approximation [Eq. ([T])] and for this reason it can 
hardly be applied to strongly interacting systems. 

In this paper we propose a more physical way of looking 
at atomic shells and bonds based on the idea of quantum 
stress focusing, by which we mean the spontaneous for- 
mation of regions of high and low pressure in the electron 
gas. The existence of such regions is nontrivial in view 
of the fact that the electrostatic potential created by the 
nuclei has no maxima or minima (Earnshaw's theorem). 

To understand how pressure maxima can arise we must 
recall a few basic concepts from the mechanics of contin- 
uum media. First, we introduce the stress tensor [l^, 11 1, 
Pij{r), - a symmetric rank-2 tensor with the property 
that pij is the i-th component of the force per unit area 
acting on an infinitesimal surface perpendicular to the 
j-th axis [l3|- The divergence of this tensor, J2j djPij{r) 
{dj being the derivative with respect to rj , the j-th com- 
ponent of r), is the i-th component of the force per unit 
volume exerted on an infinitesimal element of the elec- 
tron gas by the surrounding electrons. For convenience, 
in the following we will exclude from the stress tensor the 
electrostatic Hartree contribution, which is then treated 
as an external field on equal footing with the field of the 
nuclei. 
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FIG. 1: Indicators of electron localization as functions of 
the radial distance r constructed respectively from the non- 
interacting and interacting pressures and p [Eq. ((9)l] for 
the Ar, Zn, and Xe atoms. The XC part of the pressure is 
evaluated with the GGA functional of Ref . [l^ . The maxima 
of the pressure reveal electronic shells, indicated as usual by 
the capital letters K,L,M,... 

Next, we separate pij into a hydrostatic pressure part, 
which is isotropic, and a shear part, which is traceless: 

Pij{r) = Sijp{r) + TTijir), Tr tt^- = 0. (3) 

Here p{r) = ^Tr Pij{r) is the quantum pressure, i.e. the 
derivative of the internal energy with respect to a local 
deformation that changes the volume of a small element 
of the electron liquid without changing its shape. Sim- 
ilarly, the shear stress 7rij(r) is the derivative of the in- 
ternal energy with respect to a local deformation that 
changes the shape of a small element of the liquid, with- 
out changing its volume. In a quantum system, these lo- 
cal deformations are implemented through a local change 
in the metrics (see Eq. (ITU) below and Ref. 11 for details). 



In the ground state of an electronic system the forces aris- 
ing from the divergence of pij must exactly balance the 
external force exerted by the nuclei. In other words, the 
stresses satisfy the equilibrium condition 

- nir)dMr) = d^r) + ^ d,n,, (r), (4) 

j 

where v{r) is the external potential due to the nuclei plus 
the Hartree potential. 

Eq. ^ is not very useful if we don't have explicit ex- 
pressions for p and iTij . However, certain approxima- 
tions can be made. For example, replacing the exact 
quantum pressure by the quasi-classical Fermi pressure 
p^^{r) = |(37r^)^/'^n^/'^ and neglecting the shear term 
yields the well-known Thomas-Fermi equation for the 
equilibrium density. Because the electrostatic potential 
v{r) in an atom does not admit maxima and minima, we 
immediately see that this approximation (as any approx- 
imation that neglects n) cannot yield local maxima or 
minima in the pressure. In the special case of the spheri- 
cal Thomas-Fermi atom this implies that the density and 
the pressure are both monotonically decreasing functions 
of the radial distance r: the shell structure is absent. 

The situation changes radically when the shear stress 
is included. Now it is possible for the hydrostatic force 
dip to vanish because the electrostatic force ndiV can 
be balanced by the shear force '}2,j dji^ij ■ It is there- 
fore possible to have maxima or minima in the pressure. 
We may picture the regions of high pressure as regions 
which are hard to compress, because of the high en- 
ergy cost of bringing the particles closer together against 
the "exchange-correlation hole" . These "maximally com- 
pressed regions" quite naturally correspond to the shells 
and bonds of descriptive chemistry. 

We now show that regions of maximum and minimum 
pressure do occur in spherical atoms, with peaks corre- 
sponding to the K,L,M... shells of the standard shell 
model. To see this we don't need to go any further 
than the lowest-order approximation for the stress tensor, 
which has the form: 

P% = \ll {d,ri.d,^l>i, + d,rMi^) - i% V^n, (5) 

la 

where the sum runs over the occupied Kohn-Sham or- 
bitals. For an atom of spherical symmetry this further 
simplifies to 

P^J = P ^r] + j ' 

where p^ is the nonintcracting pressure given by 

la 

and 

^^ = Y.{\dr^ia?~r-^\de^i,?). (8) 

la 
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FIG. 2: Non-interacting pressure (p^) and shear stress (vr^) 
pressure in Ar. The peaks in identify the position of the 
shells and the peaks in tt^ define the boundary between the 
shells. Inset: distribution of stresses in the atom. The grey 
scale indicates the pressure distribution with dark (light) re- 
gions corresponding to high (low) value of the pressure in- 
dicator L (Eq. [9]). The arrows indicate the direction of the 
shear force (2/3)V7r + {2/r)n. Observe how the shells are 
"squeezed" by shear forces pointing in opposite directions. 

To graphically represent the noninteracting pressure 
, we define the dimensionless quantity — p^ /p^^ . 
Since this ratio diverges both at a nucleus (with a positive 
sign) and in the density tail (with a negative sign) [l^, 
we find it convenient to plot the function 

L{r) = -\ l+ 1 (9) 

which has values in the range (0, 1), with L = \ cor- 
responding to p^ = 00 and L = corresponding to 
p^ = — oo. Figures la-lc show conclusively the exis- 
tence of peaks and troughs of the noninteracting pres- 
sure (solid curves), each of which defines a spherical sur- 
face on which the gradient oip^ vanishes. Since p^^ is a 
monotonically decreasing function, the oscillatory behav- 
ior of these graphs is entirely due to p^: in particular, the 
maxima and minima of L are very close to the maxima 
and minima of p^ . As discussed above, we interpret the 
peaks in pressure as quantum shells. Figure 2 shows the 
behavior oip^ and tt^ as functions of r, and the inset de- 
picts schematically the distribution of shear forces in the 
atom. Observe that the minima of the shear stress are 
in correspondence with the maxima of the pressure, and 
viceversa. This is natural, because the gradients oip and 
TT, i.e., approximately, the bulk force and the shear force 
must have opposite signs and largely cancel in order to 
balance the relatively weak nuclear attraction. The fact 
that the shell is located at the minimum of the shear 
stress implies that the shear force has opposite signs on 
the two sides of the shell, thus "locking" the shell into 
position. 

Two more points should be made about Eq. ([7]) for the 
pressure. The first is that the Laplacian of the density 



appears in it prominently, as a universal component of 
the pressure: this explains a posteriori the partial suc- 
cess of the Laplacian as an indicator of shell structure. 
More interestingly, the ELF (Eq. [2]) is also closely related 
to the pressure. Indeed, the non-interacting pressure in- 
cludes a "bosonic" contribution, pq, which is obtained 
from Eq. ([7]) by putting all the electrons in the lowest 
energy orbital "tpoa'- 

poa = (l/3)[|Vn,|V(4n.) - (3/4) V^n,]. (10) 

The excess pressure pcxa = p^ — poa (r) is called "Pauli 
pressure" (since it is generated by the Pauli exclusion 
principle, which forces the occupation of higher energy 
orbitals) and is easily seen to coincide with D^, the cur- 
vature of the conditional probability function (Eq. ([1])) 
and main ingredient of the ELF. Physically, high elec- 
tronic pressure means that it is difficult to bring the par- 
ticles closer together, which implies that the XC hole 
surrounding the electron is very "deep" , i.e. the proba- 
bility of finding another electron of the same spin around 
the reference electron is low. This leads to a high value 
of the ELF. 

One of the attractive features of our proposal is the 
relative ease with which strong electron-electron inter- 
actions can be included in the calculation of the stress 
tensor. Already p^ includes some interaction effects 
through the Kohn-Sham orbitals [l3|. To include addi- 
tional exchange-correlation (XC) effects we observe that 
in density functional theory (DFT) the XC energy is a 
functional not only of the density n, but also, tacitly, of 
the metrics gij of the space[ll|: E^c = £'xc['^, 

In conventional applications of DFT, one works with 
fixed Euclidean metrics gij = Sij, so there is no need to 
emphasize the dependence on the metrics. Here, how- 
ever, we are interested in small deformations, which can 
be described by the coordinate transformation ^(r) = 
r — u(r), where u is an infinitesimal nonlinear function 
of r. The transformation from r to ^ changes the density 
from n(r) to, up to first order in u, n(^) = n{^)+Uidin{^) 
(notice that i here are cartesian indices with the conven- 
tion that repeated indices are summed over), and the 
metric tensor from the Euclidean value Sij to 

gij = {dra/d£,i){dra/d^j) = Sij + diUj + djUi. (11) 

Since the XC energy cannot depend upon the arbi- 
trary choice of coordinates we must have Ey^c[n,5ij\ — 
Eyic[n.,gij\, which for an arbitrary infinitesimal transfor- 
mation implies 

na.^ = 9,(^^-2^1, (12) 

where the derivative with respect to n is taken at con- 
stant gij and viceversa. Since the left-hand side of this 
equation is the XC force density, the right-hand side must 
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be the divergence of the XC stress tensor. This leads to 
the result 



p7 



(13) 



where Vxc is the XC potential defined by v^c = SE^c/Sn, 
and the functional derivative is evaluated at gij — Sij. 
This important result shows that we may calculate the 
XC contribution to the stress tensor from a knowledge of 
the XC energy as a functional of density and metrics. It 
turns out that semi-local XC functionals 3i 12 1 are of a 
form that can be easily generalized to include a nontrivial 
metrics. 

For example, consider the generalized gradient approx- 
imation (GGA) for a spin-unpolarized system. The stan- 
dard form of this functional is i?xc["-] = / dre^dn, s), 
where s = \\/n\/{2kFn) is the reduced density gradient. 
Going to curvilinear coordinates we immediately obtain 
the generalized form 



(14) 



where g = det(5ij) is the determinant of the metric 
tensor, n(^) = n(r(^)) is the density in the new co- 
ordinates, and s = g^^ dihdjfi/ (2/ci?n), where kp ~ 
(STT^n)^/^ and g'-' = {g~^)ij. Using the identity 5g = 
gg^-'Sgij = —ggtjSg^^ we evaluate the functional deriva- 
tive of Exc[n, gij] with respect to g^j, and then substitute 
it into Eq. (fTS]) . We finally obtain 



diudjU 



2fci?n|Vn| ds 



The trace of this yields the XC pressure 

s de^c 



P 



3 ds 



(15) 



(16) 



which exactly recovers the uniform-gas limit |10 |. 

In Figs. la-lc we have also plotted the normal- 
ized pressure L{r) including the XC correction (dashed 
curves), which is calculated by the GGA of Ref. [l6|. In 
the core region of the atom, where the density is high, p^'^ 
is utterly negligible. In the valence region, the XC con- 
tribution becomes relatively important, so the difference 
is noticeable. In the density tail, both and p^^ decay 
exponentially, causing the difference to be small. It is 
possible, however, and indeed quite likely, that XC con- 
tributions become more important in strongly correlated 
systems, where the zeroth-order description provided by 
the Kohn-Sham orbitals begins to break down. This re- 
mains a subject for future investigations. 

In conclusion, we have found that the analysis of the 
quantum stress field reveals important features of the 
electronic structure of atoms. Although we have con- 
sidered only atoms so far, there is every reason to expect 



that stress focusing will also occur in molecules in cor- 
respondence of electron pair domains, such as the ones 
that are commonly associated with single and multiple 
covalent bonds and with the so-called lone-pair regions 
of the molecule. In general, the local pressure - the trace 
of the stress tensor - is expected to provide an excellent 
quantitative description of electron localization. More 
information may emerge from the complete study of the 
stress tensor (its eigenvalues and principal axes) and its 
topology - a study that is just beginning to be under- 
taken. 
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